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Abstract. - We investigate a self-organized swimmer at low Reynolds numbers. The microscopic 
swimmer is composed of three spheres that are connected by two identical active linker arms. 
Each linker arm contains molecular motors and elastic elements and can oscillate spontaneously. 
We find that such a system immersed in a viscous fluid can self-organize into a state of directed 
swimming. The swimmer provides a simple system to study important aspects of the swimming 
of micro-organisms. 
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Introduction. — Micro-organisms often use flagella 
or cilia to move autonomously in an aqueous environ- 
' ' ment [1]. In eukaryotes these hair-like appendages are 
*J internally driven by molecular motors. The motors set 
£f-}the appendage into oscillatory motion leading to two- or 
£N| three-dimensional beating patterns, which propel the or- 
ganism. Due to the complicated structure of cilia and flag- 
ella, our understanding of swimming through these intcr- 
(~* > nally driven filaments is still far from being complete. This 
ON holds even more for effects caused by hydrodynamic in- 
^^tteractions between swimming micro-organisms, which can 
lead to turbulent motion in bacterial suspensions [2] or to 
the formation of vortex arrays for ensembles of sperma- 
^tozoa [3]. In the latter case, the flagella propelling the 
CS Spermatozoa have been observed to beat synchronously. 
Synchronous beating can lead to metachronal waves [4] 
which are of biological importance for the swimming of the 
single-celled freshwater ciliate Paramecium or the trans- 
port of mucus in our lungs. 

Physical studies of the swimming of micro-organisms 
have focused on the motion resulting from a given se- 
quence of shape changes. As the Reynolds numbers asso- 
ciated with the flow generated by such swimmers are low, 
inertial effects do not play a role. Consequently, swim- 
ming requires shape changes of more than one degree of 
freedom [5] . Different swimmer geometries have been con- 
sidered in this context, see, for example, Refs. [5-10] the 
simplest being probably the swimmer introduced by Najafi 
and Golestanian [9,11-13], which consists of three lined-up 
spheres with periodically varying distances, see Fig. la. 
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Fig. 1: a) Schematic representation of the three-sphere system, 
b) Schematic representation of the dynamic linker element. It 
consists of a filament, a motor filament, and a linear spring with 
elastic modulus K. Motors bind to the filament with rate uib 
and unbind with rate u> u . Motors bound to the filament move 
with velocity v so as to shorten the linker length x. Motors 
are connected to the backbone by springs of stiffness k and 
extension y along the z-axis. 



The mechanisms underlying these shape changes, how- 
ever, have received less attention in the physical com- 
munity. Exceptions are provided by works on the os- 
cillatory beating of flagella [14-16]. There, the periodic 
motion is assumed to be generated by spontaneous me- 
chanical oscillations of molecular motors coupled to elas- 
tic elements [17, 18]. Combining the action of motor 
molecules and the bending elasticity of a flagellum, wave 
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patterns similar to those observed experimentally arc gen- 
erated [14-16]. 

In this work, we introduce a simple dynamical system 
consisting of three spheres that are arranged in a line 
and that are connected by linkers. The linkers consist 
of molecular motors and clastic elements and their length 
can oscillate spontaneously. While the swimmer is sym- 
metric with respect to space-inversion the system can dy- 
namically break this symmetry and spontaneously swim 
dircctionally. We investigate the system's bifurcation di- 
agram and provide evidence that the transition to direc- 
tional swimming shares striking similarities to transitions 
occuring in integrate-and-fire neurons. 

The active three-sphere swimmer. Consider 
three spheres, each of radius R moving in a fluid of vis- 
cosity 7/ along the x-axis, with neighboring spheres being 
linked by elements as illustrated in Fig. lb. Their struc- 
ture resembles that of sarcomeres, which are the elemen- 
tary force generating units of skeletal muscle, and consist 
of motors moving along a polar filament combined with 
a linear spring. Such structures are known to be able to 
oscillate spontaneously [18], a phenomenon that has been 
observed experimentally for muscle sarcomeres [19]. Be- 
fore investigating the dynamics of the three actively linked 
spheres, let us first describe the three-sphere swimmer 
with externally imposed changed of the spheres' mutual 
distances. 

The three spheres are located at positions Xj, i — 1, 2, 3 
along the x-axis. At low Reynolds numbers, the forces fj 
acting on sphere j with j = 1,2,3 result in velocities Xj 
according to [20] 

x i = ^H ij f j . (1) 

j 

Here, the mobility tensor H depends on the fluid's viscos- 
ity, the size of the spheres, as well as the distances between 
them, and satisfies ££y = Hji. Using Eq. (1) as well as 
the global force balance equation fa + fa + fa = 0, the ve- 
locity of the center of mass of the system, dx s /dt = x s = 
(xi + X2 + X3)/3, can be expressed in terms of the dis- 
tances between the middle and the outer spheres, x\ and 
x T , respectively. For periodically varying distances, one 
finds for the displacement x s after one period T, 

x s = x s dt = C{x\, x T ) dxi A dx r , (2) 
Jo Jo 

where O denotes the oriented surface encircled in the 
(sci, x r )-plane during one period of the motion, see Fig. 2a 
and also Refs. [5,7]. In this expression, C is a field with 
C(x\,x T ) = C(x T ,x\) that only depends on the mobility 
tensor H. Thus the displacement is independent of the 
oscillation period. For the calculation of C, it is conve- 
nient to employ the Oseen approximation, which consists 
in expanding H in terms of R/x\ and R/x r and retaining 
leading order terms only. In lowest order, this corresponds 
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Fig. 2: a) Schematic representation of possible orbits in the 
(x'i, Xr)-plane of periodic motions of the three-sphere system. 
Left: reciprocal motion, middle: non-reciprocal motion with 
x s = 0, right: non-reciprocal motion with x s 7^ 0. b), c) 
Left and right linker length as a function of time in the non- 
swimming state Oi, see Fig. 4, with £ > £ s and d), e) in the 
swimming state Sw with £ < £ s . Other parameters as in Fig. 3. 

to neglecting hydrodynamic interactions between different 
spheres. In this case, C(xi,x r ) = 0, implying that hy- 
drodynamic interactions are essential for swimming of the 
three-sphere system. 

If motion is reciprocal, that is, if the corresponding 
curve in the (x\, x r )-plane retraces itself, then x s = illus- 
trating Purcell's Scallop theorem, see Fig. 2a left. Note, 
that a trajectory which is invariant under the exchange of 
left and right, i.e., space inversion, is necessarily recipro- 
cal. Examples of other periodic motions with vanishing 
net displacement are given by x\(t) = x T (t + T/2), see 
Fig. 2a middle. These solutions are invariant under space 
inversion and a simultaneous shift in time by T/2. A sim- 
ple form of x\ and x r with x s ^ was proposed by Najafi 
and Golestanian [9], see Fig. 2a right. 

We now describe the dynamics of the linker elements de- 
picted in Fig. lb, consisting of motors and elastic compo- 
nents. Similar systems have been studied in the contexts 
of spontaneous spindle [21] and muscle oscillations [22]. 

A linker element consists of motors moving dircctionally 
along a polar filament and a linear spring of stiffness K. 
The motors are attached to a common backbone by springs 
of stiffness k and extension y. In a sarcomere, the motors 
would be ensembles of myosin II that move on an actin 
filament, while elastic elements are provided by various 
structural elements, e.g., the protein titin. Motors bound 
to the filament move such that the linker shortens. An un- 
bound motor binds to the filament at rate wt>- Attached 
motors unbind from the filament at rate u) u . These rates 
generally depend on the force applied to the motor. We 
assume that the force dependence is restricted to the un- 
binding rate. Motivated by Kramers' rate theory we write 
uj u = exp{|/| a/fceT}, where f = ky and a is a mi- 
croscopic length scale. Also the velocity of bound motors 
depends on the applied force. For simplicity, we assume a 
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linear force- velocity relation v(f) = Vq (1 — f/fa)- Here, fa 
is the stall force and vq the velocity of an unloaded motor. 

The dynamic equation for the length x of the linker 
element connecting two spheres, say 1 and 2, x — x 2 —x\, is 
given by Eq. (1), where the force fa on sphere j is the sum 
of the elastic and the motor forces generated by the linker. 
The elastic force fa is given by fa = —K{x — Lq) with 
I/O being the spring's rest length. To calculate the force 
exerted by the motors, / m , consider first a single motor 
bound to the polar filament. If y n denotes the extension 
of the spring linking this motor to the backbone then the 
force it exerts is —ky n . The value of y n changes according 
to 

Vn = V n + X , (3) 

where v n is the velocity of this motor on the filament. The 
total motor force then is f m = —k^2 n=1 (T n y n . Here, N 
is the overall number of motors and a n = 1 if motor n is 
bound and otherwise. We employ a mean-field approxi- 
mation, which consists in assuming that all bound motors 
have the same spring extension, y n = y for all n. We de- 
termine the mean position by the product of the average 
motor speed multiplied with the average time they stay 
attached, y = y/u> u . For unbound motors, we assume fast 
relaxation, such that their distribution equals the equilib- 
rium distribution implying y = for these motors. Under 
these conditions, the fraction Q of bound motors evolves 
according to [21] 

Q = - (wb + w u ) Q , (4) 

while the total motor force is 

fm = -N(x)Qky . (5) 

Here, N(x) — (£f + £ m — x)/A is the number of motors in 
the overlap region of the filament and the motor backbone 
of lengths If and £ m , respectively, and A is the distance 
on the backbone between adjacent motors. 

Before considering the active three-sphere swimmer, let 
us present the dynamics of two actively linked spheres. 
The dynamic equations for this system are given by 
Eqs. (1), (3)- (5), where the forces on the spheres 1 and 
2 are fa = —fa = —fa — fm- To lowest order the Oseen- 
expansion of H gives Hn = H 22 = (6 7r 77 1 , and 
H 12 = 0. Ha is the Stokes friction of an isolated sphere. 
In dimcnsionlcss form, the equations of motion read 

i'x = -2- [N(x)Qny + x-L] = -2F(x,Q,y) (6) 
Q = l-[l + w(y)]-Q (7) 
x = y-K£) + 7 ]-l , (8) 

where the dimensionless parameters are given by £ = 
6Trr)Ru)b/K, L = Lqlo^/vq, k = k/K, ui = uJu/lo^ and 7 = 
kvo j (fo^b) ■ Furthermore, x = xu>i,/vq and y = yuj-^/vo, 
while time has been rescaled by ujb- 

The equations have a stationary state (Jo, Qo, yo)- For 
sufficiently short times, a small perturbation 5x, 5Q, Sy 



will evolve as Sx(t) oc exp(si) and analogously for 6Q 
and Sy. It grows if 5R(s) > 0. This condition can be 
expressed in terms of the system parameters 1 and implies 
yau)' u (yo) > lo^ + <Ai(yo)- ^ ne f° rce dependence of the 
binding kinetics is thus essential for an instability of the 
stationary state. The system oscillates if 3(s) 7^ at 
the instability. This implies Qq kyo < K A, that is, the 
increase of the motor force by adding one motor in the 
overlap region needs to be smaller than the corresponding 
increase of the restoring force by the elastic element to 
produce an oscillatory instability. 

The oscillation mechanism can be understood intu- 
itively as follows: As the motors shorten the linker, the 
elastic restoring forces increases and, therefore, the force 
on the motors. This in turn increases the unbinding rate of 
motors. When a few motors detach from the filament, the 
remaining motors experience an even higher force, such 
that an avalanche of motor unbinding events occurs. The 
elastic element then stretches the linker and the motors 
rcbind, repeating the cycle. 

This behaviour shares similarities with the dynamics of 
integrate-and-fire model neurons [23]. There, the electric 
membrane potential of a neuron increases monotonically 
as a consequence of an input current. Upon reaching a 
threshold value, the potential is instantaneously reset and 
the loading process restarts. In the case of the linker ele- 
ment, the mechanical analog of neuron's electric potential 
is the internal stress of the linker. Loading occurs through 
the action of motors and resetting is accomplished by the 
spring K after a detachment avalanche of the motors. 

Let us now return to the three-sphere system, see Fig. 1, 
where the spheres are connected by dynamic linker ele- 
ments. The forces acting on the three spheres are /1 = 
-fe - /m> h = ~h ~ /a, and fa = f% + where the 
superscripts I and r distinguish between the left and the 
right linker, respectively. 

As we have seen above, hydrodynamic interactions be- 
tween different spheres are essential for swimming of a 
three-sphere system. Compared to the two-sphere case, 
we thus now expand the mobility tensor to the next higher 
order. This leads to Hn = (6 7r 77 — 1 , i = 1,2,3, 
H 12 = (4 7rryi?)- 1 i?/x 1 , H 23 = (4tt?7 R)- x R/x t , and 
H13 = (i'7T7'iR)~ 1 R/(xi + x r ). The equations of motion 
are then given by Eqs. (7) and (8) respectively for the left 
and right part of the swimmer. Equation (6) is replaced 

by 

££i,r = {¥- -2\F(x hl ,Qi tI ,y l , I ) (9) 

[2 X\ Xr J J 

with r = Ruj^/vq. As before, the dynamic equations are 
invariant under space inversion, i.e., exchange of left and 

lr The instability occurs at f(l + u(y)) + 1 - QKyA+ [N(x)Qk(1 + 
v(y) - yu'(y)]/[w(y) + yw'(y) + 7] = 0, with A = v /Aw b . Here, x, 
Q and y are taken in the stationary state. 
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Fig. 3: State diagram of the dynamic three-sphere swimmer. 
For £ > £1 the system is stationary. Filled diamonds indicate 
the mean velocity of the swimmer's center of mass. For £ s < 
£ < £1 it oscillates, but does not move on average, v a = 0. For 
£ < £ s , v s > 0. Crosses indicate the relative phase between the 
oscillations of the two linker elements. Equations (7)-(9) were 
solved numerically with L = 137.5, n = 400, u>o = 0.5, 7 = 9.1, 
A = 0.16 and r = 11. 



Fig. 4: Bifurcation scenario represented by the linker ampli- 
tude i max for varying parameter £. Solid and dotted lines 
indicate stable, dashed lines unstable states. For £ < £1 the 
stationary state St becomes unstable via a Hopf-bifurcation. 
The emerging oscillatory state Oi loses stability at £ = £ s and 
the swimming solution Sw (dotted line) appears. Furthermore, 
at £ = £2, a second oscillatory Hopf-mode O2 appears. Inset: 
Representation of the swimming bifurcation in terms of the 
I/2-norm (see text). 



right. 

Figure 3 summarizes the three-sphere system's behav- 
ior with the active linker elements as the parameter £ is 
changed. This can be achieved, for example, by chang- 
ing the viscosity of the surrounding fluid or the binding 
rate of the motors. Having £ > £1, the system settles 
into a stationary state. For £ < £1, the linkers sponta- 
neously oscillate in length. The oscillations of the two 
linker elements arc identical but shifted with respect to 
each other by half a period. As is illustrated in Fig. 2b, c, 
this symmetry implies that, for low Reynolds numbers, 
the swimmer's center of mass x s remains stationary on 
average. There is a second critical value £ s , at which the 
relative phase shift between the two oscillating linkers <f> 
starts to deviate from half a period, such that the system 
spontaneously swims if £ < £ s . In contrast, the form of 
the oscillatory motion of each of the two linker elements 
essentially does not change. By dynamically breaking the 
space-inversion symmetry, the corresponding orbit in the 
(xi, x r )-planc is not symmetric with respect to the line 
X \ — Xj^ } SGG Fig. 2d,e, hydrodynamic interactions allow in 
average for directed motion of the three-sphere system. 

Decreasing the value of £ further speeds up the swim- 
mer. The increase in the swimming velocity is the conse- 
quence of two effects. On one hand, the relative phase shift 
gets larger and, on the other hand, the amplitude of the 
oscillation increases, both leading to larger areas encircled 
by the orbit in the (aq, £ r )-planc. For still lower viscosi- 
ties one expects eventually a decrease in the velocity such 
that swimming stops at £ = 0. The dynamic equations 
presented here do not capture this effect as their applica- 
bility is restricted to the regime of low Reynolds numbers. 

By the symmetry of the swimmer, with each moving 



state coexists a state moving into the opposite direction. 
The actual direction of motion is determined by the initial 
conditions. Let us also note that changing the parameters 
7 or r can also lead to transitions from the stationary 
to the non-swimming oscillatory state or from the non- 
swimming oscillatory to the swimming state. We did not 
find two subsequent bifurcations leading from the station- 
ary to the swimming state, though. 

We will now take a closer look at the swimming tran- 
sition. The bifurcation diagram of the swimmer in terms 
of the oscillation amplitude i max is presented in Fig. 4. 
The stationary state St loses stability through a super- 
critical Hopf-bifurcation at £ = £1. The newly gener- 
ated oscillatory state Oi satisfies x\(t) = x r (t + T/2), 
where T is the dimcnsionlcss oscillation period. At £ = 
£ s , Oi loses stability through a pitchfork bifurcation, 
see Fig. 3, resulting in the swimming state Sw. Here, 
x\(t) = x r (t + </>), where the relative phase </> ^ T/2. For 
£ < £ S ; the states Oi and Sw have almost the same am- 
plitude. Therefore, we display in the inset of Fig. 4 a 
representation of the bifurcation in terms of the L 2 _n orm, 
T- 1 J T dt (xf + xl + Qf + Q 2 r + yf + . For a similar 
value of £ a second non-swimming state O2 bifurcates 
from the stationary state at £ = This state, which 
is unstable for all values of £, satisfies x\(i) = x T (t). 

The appearance of the state O2 together with an in- 
spection of the corresponding isotropy lattice might sug- 
gest mode interaction to be at the origin of the symmetry 
transition. However, the frequencies of and 2 are 
different at £ s , which argues strongly against mode inter- 
action. 

While we are currently lacking a thorough understand- 
ing of the swimming transition, we would like to point out 
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that this transition is reminiscent of a phenomenon occur- 
ing for two coupled integrate-and-fire neurons [24]. The 
membrane potential of an integrate-and-fire neuron is in- 
creased at a constant rate and resetted when a certain 
threshold voltage is reached. Upon resetting, the neuron 
sends an electric pulse to other neurons it is connected 
to and influences their membrane potentials. For broad 
pulses, the two coupled neurons settle into a state in which 
they fire periodically with a relative phase shift of half the 
period. For short enough pulses, however, this state is un- 
stable and a new value of the phase-shift is generated [24] . 
For the swimmer, in analogy to the electrical pulses for 
the neurons, mechanical pulses are transmitted netween 
the two linkers. A mechanical pulse is released by one 
linker when it relaxes as a consequence of a detachment 
avalanche of the motors. In that case, the duration of the 
pulse changes, because of alterations of the duration of the 
linkers expansion phase that effectively transmitts the me- 
chanical pulse. Decreasing £ leads to faster stretching of 
the linkers. In terms of physical parameters, decreasing £ 
either corresponds to reducing friction forces on the beads, 
to decreasing the binding rate of motors, or to strength- 
ening the elastic elements. As for the two-neuron system, 
a sharper pulse can be sufficient to cause dephasing [25], 
which implies swimming. 

The equations of motion (7)-(9) are of the form of a 
singularly perturbed system with £ as the small parame- 
ter. Such systems are known to present strong deforma- 
tions of a limit cycle close to a Hopf-bifurcation, where 
the deformation occurs for exponentially small parameter 
changes [26]. In certain parameter regions, we indeed find 
that the Hopf-bifurcation is followed by such a deforma- 
tion which is called a Canard explosion. Accompanying 
with the Canard phenomenon, we find also excitable dy- 
namics for the swimmer. 

In addition to the Canard explosion, the singularly 
perturbed system (7)- (9) shows dynamic behavior that is 
reminiscent of the dynamics associated with global bifur- 
cations. In this case, the spontaneous oscillations of the 
linker are still generated by a Hopf bifurcation. The dy- 
namics close to the bifurcation, however, shows character- 
istics of the dynamics close to a saddle-node infinite period 
bifurcation (SNIPER). A SNIPER bifurcation is charac- 
terized by two fixed points, a saddle and a node, that 
collide, which leads to oscillations around a third, unaf- 
fected fixed point, namely, an unstable node [27]. In the 
present case, this limit cycle is symmetric with respect to 
space-inversion. When decreasing the control parameter £ 
further, the limit cycle splits up into two limit cycles that 
are mirror images of each other upon space inversion. This 
transition is similar to an anti-gluing bifurcation [28]. A 
detailed analysis of these bifurcations will be given else- 
where. 

Discussion. — In this work, we have presented a sim- 
ple self-organized swimmer. It consists of three aligned 
spheres connected by active linker elements that can os- 



cillate spontaneously. By decreasing the dimcnsionlcss pa- 
rameter £ below a critical value, which can be achieved ei- 
ther by reducing the the viscosity of the surrounding liquid 
or the motors' binding rate or by increasing the stiffness K 
of the elastic clement, the system oscillates spontaneously. 
While this state is not associated with an average displace- 
ment, a subsequent bifurcation encountered when decreas- 
ing £ further, leads to directed swimming. Changing other 
system parameters can either lead to spontaneous oscilla- 
tions or to swimming but has not been found to induce 
both transitions subsequently. 

Could the self-organized three-sphere be realized exper- 
imentally? One approach could be based on muscle sar- 
comeres as linker elements. As mentioned above, sarcom- 
eres can oscillate spontaneously [19]. Using parameter val- 
ues appropriate for sarcomeres 2 and assuming that they 
link spheres with a diameter of R = 0.1 /zm, we find swim- 
ming velocities on the order of v s w 1 /im/min. This has 
to be compared to the diffusion constant. For a sphere 
with a diameter of 0.1/xm we find in water at room tem- 
perature D w 2/im 2 /s. On a distance comparable to the 
swimmer size L w 3 fim, diffusion thus dominates the 
motion Lv s /D = 0.025. In order to reach bigger swim- 
ming speeds, the oscillation frequency could be increased. 
To this end stiffer springs and stronger motors would be 
necessary. Another way to increase the swimming speed 
would be through modifications of the swimmer design. 
For example, it has been argued that the pushmepullyou 
swimmer [10] is faster than the three-sphere swimmer. 

Our self-organized swimmer can be used to study the ef- 
fects of hydrodynamic coupling between micro-swimmers. 
So far, studies of this subject have largely neglected the 
mechanism used by micro-swimmers to generate the shape 
changes necessary for swimming. Preliminary results sug- 
gest interesting synchronization effects that greatly alter 
the system's behavior compared to swimmers with fixed 
sequences of shape changes. Obviously, it would also be 
interesting to extend our approach to study the dynamics 
of cilia and flagella. 

* * * 

We thank A. Hilfinger and E. M. Nicola for useful dis- 
cussions. 

REFERENCES 

[1] D. Bray, Cell Movements: From Molecules to Motil- 
ity, Garland, New York, 2nd edition (2001). 

[2] C. Dombrowski et al, Phys. Rev. Lett. 93, 098103 
(2004). 

[3] I.H. Riedel, K. Kruse, and J. Howard, Sci- 
ence 309, 300 (2005). 



2 i { = 0.6 /Ltm, £ m = 0.75 Aim, A = 56 nm, k = 4pN/nm, 
K = O.OlpN/nm, L = 1.25 /jm, v = 2/jm/s, f = 4pN, 
k B T = 4pNnm, a = 3nm, uj h = 220/s and uj° = 120/s. 



S. Giintherl,2 K. Krusel,2 



[4] M.A. Sleigh, in Cilia and Flagella edited by M.A. 

Sleigh, (Academic Press, London) 1974, p. 287 
[5] E.M. Purcell, Am. J. Phys. 45, 3 (1977). 
[6] G.I. Taylor, Proc. R. Soc. London Ser. A 209, 447 

(1951). 

[7] A. Shapere and F. Wilczek, Phys. Rev. Lett. 58, 
2051 (1987). 

[8] E.M. Purcell, Proc. Natl. Acad. Sci. USA 94, 11307 
(1997). 

[9] A. Najafi and R. Golestanian, Phys. Rev. E 69, 

062901 (2004). 
[10] J.E. Avron, O. Kenneth, and D.H. Oaknin, New 

J. Phys. 7, 234 (2005). 
[11] R. Golestanian and A. Ajdari, Phys. Rev. E 77, 

036308 (2008) 
[12] R. Golestanian, Eur. Phys. J. E 25, 1 (2008) 
[13] R. Golestanian and A. Ajdari, Phys. Rev. 

Lett. 100, 038101 (2008) 
[14] S. Camalet, F. Julicher, and J. Prost, Phys. Rev. 

Lett. 82, 1590 (1999). 
[15] I.H. Riedel-Kruse, A. Hilfinger, J. Howard, and 

F. Julicher, HFSP Journal 3, 192 (2007) 
[16] S. Camalet and F. Julicher, New J. Phys. 2, 24 

(2000). 

[17] C.J. Brokaw, Proc. Natl. Acad. Sci. USA 72, 3102 
(1975). 

[18] F. Julicher and J. Prost, Phys. Rev. Lett. 78, 4510 
(1997). 

[19] N. Okamura and S. Ishiwata, J. Muscle Res. Cell 

Motil. 9, 111 (1988). 
[20] G.K. Batchelor, J. Fluid. Mech. 79, 1 (1976). 
[21] S. Grill, K. Kruse, and F. Julicher, Phys. Rev. 

Lett. 94, 108104 (2005). 
[22] S. Gunther and K. Kruse, New J. Phys. 9, 417 

(2007). 

[23] L. Lapicque, J. Physiol. Pathol. Gen. 9, 620 (1907) 
[24] C. van Vreeswijk, L.F. Abbott, and G.B. Er- 

mentrout, J. Comput. Neurosci. 1, 313 (1994) 
[25] M. Timme, private communication. 
[26] E. Benoit, J.L. Callot, F. Diener, and M. Di- 

ener, Collect. Math. 32, 37 (1981) 
[27] S.H. Strogatz, Nonlinear Dynamics and Chaos, 

Addison- Wesley, (1994) 
[28] P. Coullet, J.M. Gambaudo, and C. Tresser, CR 

Acad. Sci. I-Math. 299, 253 (1984) 



